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Advanced structural reliability methods are utilized on the PSAM project 
to provide a tool for analysis and design of space propulsion system 
hardware. The role of the effort at the University of Arizona is to provide 
reliability technology support to this project. Reliability refers, in 
general, to application of probabilistic and statistical methods to manage 
uncertainties in engineering design. 

PSAM computer programs will provide a design tool for analyzing uncer- 
tainty associated with thermal and mechanical loading, material behavior, 
geometry, and the analysis methods used. Specifically, reliability methods 
are employed (a) to perform sensitivity analyses (a study of the importance of 
each design factor), (b) to establish the distribution of a critical response 
variable, (e.g., stress, deflection), (c) to perform reliability assessment, 
and (d) ultimately to produce a design which will minimize cost and/or weight 
(Fig. 1). 

Uncertainties in the design factors of space propulsion hardware, as 
defined in Fig. 1, are described by probability models constructed using 
statistical analysis of data. Two examples are provided in Fig. 2 and Fig. 

3. Turbine blade twist (a geometric variable) has a statistical distribution 
(Fig. 2), and significant variability in material properties are observed in 
laboratory tests (Fig. 3). Statistical methods are employed to produce a 
probability model, i.e., a statistical synthesis or summary of each design 
variable (e.g., each in Fig. 4) in a format suitable for reliability 
analysis and ultimately, design decisions. 

Given the probability model for each random variable, and the design 
equation, i.e., how the variables are related (e.g., through the NESSUS/FEM 
computer code), two basic methods for computing the distribution of a response 
variable are available. They are, fast probability integration (FPI ) , (Fig. 

4) and Monte Carlo (Fig. 5). FPI is a complex numerical algorithm which 
provides "fast" probability calculations, accurate to about ± 5% (See Refs. 1 
and 2). 

The NESSUS generated perturbed data set, an example of which is shown in 
Fig. 5, is used as a function of the input variables. Monte Carlo produces an 
exact solution, but it is expensive to run. There are a large number of 
different Monte Carlo schemes available; methods of improving efficiency of 
Monte Carlo, tailored for PSAM application, are being studied. Availability 
of the response function shown in Fig. 5 is a requirement for performing a 
fast Monte Carlo analysis. The NESSUS-generated, approximate response 
function can be used to replace the "exact" simulation models needed by Monte 
Carlo methods. 
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The result of PSAM analysis is a distribution of a response variable. But 
for design purposes, it is often important to be able to specify also the 
confidence associated with a probability statement. Confidence Intervals 
provide a quantitative description of the quality of the results of the 
probabilistic analyses (Fig. 6). Efficient methods for constructing confi- 
dence intervals are being studied. 

In summary, the University of Arizona effort on PSAM is to provide 
technical support in statistical characterization of design variables, 
advanced reliability analysis methods, and confidence interval construction 
(Fig. 7). 
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RELIABILITY METHODS ARE EFFECTIVE 
TOOLS FOR DESIGNERS 


SOURCES OF UNCERTAINTY: 

• LOADING 

• MATERIAL BEHAVIOR 

• GEOMETRY 

• STRUCTURAL MODELING 

• ANALYSIS METHOD 

• SMALL SAMPLE SIZE 



DESIGN APPLICATIONS 

• SENSITIVITY ANALYSIS 

• VARIABILITY ASSESSMENT 

• RELIABILITY ESTIMATION 

• COST/WEIGHT REDUCTION 

• CONFIDENCE LEVELS 


FIGURE 1. 
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STATISTICAL METHODS TRANSLATE 
DATA IN MODELS 



FIGURE 2. 


MATERIAL PROPERTIES EXHIBIT 
SIGNIFICANT VARIABILITY 



FIGURE 3. 
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FPI ALGORITHM IMPROVEMENTS 
DEMONSTRATED 


FATIGUE DAMAG E PR EDICTION MODEL 


• ENHANCED DISTRIBUTION LIBRARY 

- MATERIAL BEHAVIOR 

- LOADING UNCERTAINTY 

• IMPROVED ACCURACY 

• VALIDATION PROBLEMS 

- HIGHLY NONLINEAR 

- MANY VARIABLES 

• ACCURACY OF FPI i 5% 


[ X 2 1 X 2 

D 'X, + N 0 [X3<X 4 X5> 1,+ X 6 (X 4 X 5 )” 

EACH X; HAS DIFFERENT DISTRIBUTION 



• RELATIVE COMPUTER TIME 
FPI = 1 

MONTE CARLO = 60 
(CONVENTIONAL) 


FIGURE 4. 


FAST MONTE CARLO ALGORITHM (S) 
ESTABLISHED 


STRESS AT 
CRITICAL POINT 


• AN INDEPENDENT CHECK ON FPI 

• CONVERGES TO EXACT SOLUTION 



• MORE EXPENSIVE THAN FPI 

• USES NESSUS DATA BASE 


GEOMETRY 

VARIABLE 


• SEVERAL MONTE CARLO METHODS AVAILABLE 

• RELATIVE MERITS OF EACH METHOD UNDER 
STUDY 


MATERIAL 

BEHAVIOR 

VARIABLE 


NESSUS GENERATED DATA 

• DEFINES RESPONSE SURFACE 

• SUPPORTS MONTE CARLO SIMULATION 


FIGURE 5. 
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WHAT IS OUR CONFIDENCE IN THE 
ANSWER? 


PROBLEM: 


CONFIDENCE 


[■NATURAL FREQUENCY „ FREQUENCY ] , 

[of blade, n v limit. n L j 


• DIFFICULT ANALYSIS PROBLEM 

• APPROXIMATE METHOD DEFINED 



• USES NESSUS DATA BASE 

• GIVES INTERVAL ESTIMATES 

- PROBABILITY 

- VARIABILITY 


UA PROVIDES RELIABILITY 
TECHNOLOGY SUPPORT 


• LITERATURE REVIEW 

• STATISTICAL METHODS TO PRODUCE 
MODELS 

- GEOMETRY 

- MATERIAL BEHAVIOR 

• DISTRIBUTIONAL CHOICES FOR DESIGN 
VARIABLES AND RESPONSES 

• RELIABILITY ASSESSMENT METHODS 

- FPI 

- MONTE CARLO 

• CONFIDENCE INTERVALS ON 
PROBABIILTY ESTIMATES 


FIGURE 7. 
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